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ELECTRICAL RESPONSES TO TWO CLICKS: A SIMPLE 
STATISTICAL INTERPRETATION* 


/ WILLIAM J. MCGILL AND WALTER A. ROSENBLITH 


PsyCHO-AcousTIC LABORATORY, HARVARD UNIVERSITY, 
CAMBRIDGE, MASSACHUSETTS 


A probability model for the amplitude of the first neural compo- 
nent, N,, of the click response has been developed on the assumption 
that N, is the summated effect of practically simultaneous firing from 
large numbers of independent neural units. Equations derived from the 
model predict the course of recovery for the response to the second of 
two clicks. A mathematical relation is shown between the two-click sit- 
uation and the intensity growth curve for the first neural component 
of the response to a single click. This permits indirect reconstruction 
of the function relating amplitude of N, to stimulus intensity. Experi- 
ments with pairs of clicks indicate that a normal probability integral 
provides a satisfactory fit. Thus the entire intensity function of N, is 
capable of specification in terms of two statistical parameters, the mean 
and the standard deviation of the normal integral. Factors that affect 
the click response are presumed to act upon these parameters. 


When a sharp click is delivered to the ear of an anesthetized 
cat, a characteristic electrical response can be recorded from an 
electrode in contact with the round window. The response consists 
chiefly of two kinds of activity: (1) a so-called microphonic poten- 
tial, whose initial deflection is positive or negative, depending upon 
the polarity of the click, and (2) neural potentials, which are always 
negative. There are several ways of-differentiating the microphonic 
from the neural portions of the response. For our purposes, the most 
significant of these is the fact that over a wide range of intensities 
the amplitude of the microphonic is linear with the stimulus. Neu- 
ral components, on the other hand, behave in a rather complex fash- 
ion with changes in stimulus intensity. 

*This research was carried out under Contract Nb5ori-76 between Harvard 
University and the Office of Naval Research, U. S. Navy (Project Nr142-201, Re- 


port PNR-107). Reproduction for any purpose by the U. S. Government is per- 
mitted. 
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Of particular interest in connection with this paper is the earli- 
est neural component, which we shall call Ni. It is the first and 
largest neural potential recordable at the round window with a peak 
of excitation occurring approximately 1 millisecond after the acous- 
tic stimulus reaches the eardrum. Presumably it represents the re- 
sponse of neural structures near the auditory end organs (Davis et 
al., 1950). Whatever its origin, this earliest neural event is still close 
to the measurable physical stimulus and probably constitutes the 
primary neural input to the auditory nervous system. Its behavior is 
clearly related to such well-known auditory phenomena as masking 
and fatigue (Rosenblith, 1950). Thus considerable interest can be 
invested in studying the growth of the first neural component as a 
function of stimulus magnitude. 


Direct measurements of N,, however, are somewhat trouble- 
some, for the microphonic and neural components are separated in 
time for weak clicks only. They exhibit severe interaction as the 
click stimulus gets louder. Amplitude measurements of components 
of the round window response are thus likely to be considerably 
biased. At the same time experimental methods for manipulating 
the microphonic independently of the neural portions of the re- 
sponse are still in an incomplete stage of development (Kahana et 
al., 1950). In the face of such difficulties, we have attempted to re- 
construct the function relating amplitude of N, to stimulus inten- 
sity through the application of a simple probability model. 


To develop the concepts underlying the model, the click response 
must be examined in a somewhat different situation. If two clicks 
are presented to the ear, separated by a time interval of about 100 
milliseconds or less, the neural response to the second click is de- 
pressed, i.e., the ratio of the neural amplitude of the second click 
response to the resting (normal) amplitude of the same response is 
less than 1, depending upon the time separation of the clicks.* Our 
experiments indicate that the amount of depression is governed by 
three obvious factors: the intensity of the first click, the time in- 
terval between clicks, and the intensity of the second click. 


*The microphonic component of the response to the second click remains 
unaffected over the intensity range we have used. If the two clicks are separated 
by time intervals of less than 2 milliseconds, the two responses cannot easily be 
resolved and measurements of microphonics and neurals may become “con- 
founded.” For the time course of the recovery functions for the amplitude of 


the neural response to the second click see M. R. Rosenzweig and W. < 
blith (1950). g W. A. Rosen 
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We have observed too that the relative amplitude of N, in the 
response to the second click is not zero even at 2 milliseconds.* Fi- 
nally our results indicate that component N, of the second click does 
not become supernormal. 


We assume that the first neural part of the click response is 
the summated effect of firing from large numbers of neural units. 
If all neural units contribute nearly the same average voltage to the 
total response, then 

Ay =D Ri=nk, (1) 


where A vy, is the peak-to-peak amplitude of N,, R; is the voltage con- 
tributed by the ith neural unit; n is the number of neural units 


responding simultaneously, and R is the average voltage (at the 
location of the electrode). Thus it follows that the amplitude of the 
click response at the round window is proportional to the number 
of units responding. This assumption we shall discuss shortly. 


Suppose that, corresponding to any given stimulus intensity, 
there is a fixed probability that a neural unit will respond in nor- 
mal conditions. To facilitate discussion we shall label this and cer- 
tain other probabilities as follows: 


p(I) = the (resting) probability of response of a neural unit to 
the first click which is at intensity 7. 

(J) =the (resting) probability of response to the second click 
which is at intensity J. 

p*(t) = a probability operator representing the recovery of a 
neural unit that has fired. It is a function of the time since 
firing and of the intensity of the second click. 

P,. = the probability that a neural unit fires twice to successive 
click stimuli. 

Poo = the probability that a neural unit fails to fire to the first 
click but fires in response to the second. 


All possible outcomes in the two-click situation may then be repre- 
sented as follows: 


than 2 milli- 


*Listening to two clicks separated by time intervals of less selina, 


seconds confirms this statement to the extent that two clicks are p 
“more” (i.e. louder, fuzzier, more spread out) than one. 
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“Ps =p(l) p(t) pV), (2) 
Po=[1—p()]-p), (3) 
Py =p (1) -[1—p*(t) -pV)], (4) 
Po=([1—p()]1—p)I. (5) 


‘The probability, P., that a neural unit fires in response to the sec- 
ond click is the sum of Py. and Po: above: 


P,=p(1) -p* (t) -p(J) + [191 - ot). (6) 


If we postulate that p*(t) approaches zero as t, the time separation 
of the two clicks, approaches zero (i.e., in the absolute refractory 
period of any neural unit having fired in response to the first click), 
then 


fa [LS pUyien 


The relative amplitude of N, in the second click response is not zero 
unless the first click is of maximum intensity. Also we note that 
when p*(t) reaches 1 at some time t’, dependent upon the intensity 
of the second click, then: 

P,=p(J). 

t>t’ 


Let us define RA as the relative amplitude of N,, i.e., as the ratio 
of the amplitude of N, in the second click response to its resting am- 
plitude. Thus we have on the average: 


N:P. 


SS, 4 
N-p(J) ey 


where N is the total number of neural units capable of responding. 
By substituting (6) into (7) we find 


RA —1 —p(ye (ot) 1s (8) 


If we fix the intensity of the second click and the time separation 
between clicks, it is clear that RA is simply a linear function of the 
intensity of the first click. 


The experiment outlined above, in which a conditioning click 
precedes a weak test click, permits us to reconstruct the curve re- 
lating amplitude of N, to stimulus intensity. Figure 1 presents the 
actual results of such an experiment. The experimental points are 
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Figure 1. Effect of preceding click (t = 9.1 milliseconds) upon the ampli- 
tude of N, in the response to the second click. The abscissa gives the intensity 
of the first click. The second click is constant at —20 db. The curve is fitted 
by RA = 1 — ec. p(Z) in which constants are chosen as follows: ¢ = .60, 
m == —25 db, ¢ = 20 db. 


mean values of not less than 15 and not more than 25 observations 
of individual click responses from a single cat. To fit a theoretical 
curve to the data of Figure 1 we have attempted to analyze p(J) in 
a familiar way. It is well known that thresholds of individual neu- 
rons are not constant. There is indeed reasonable ground for hold- 
ing that the determinants of an instantaneous threshold are suffi- 
ciently complex to constitute the observed threshold of a neuron as 
a normal variable. The work of C. R. Pecher (1939) on sciatic nerve 
fibers in the frog confirms this and indicates that thresholds of 
adjacent neurons are independent. It is surely no accident that the 
response amplitude of the whole sciatic nerve as a function of stim- 
ulus intensity can be represented quite well by a normal probability 
integral (von Briicke et al., 1941). W. J. Crozier (1940) has argued 
at some length that variation of neural thresholds is of primary im- 
portance for understanding visual intensity phenomena. Indeed it 
was Crozier who pointed out explicitly that relations between stim- 
ulus intensity and neural response attributable to variation of neu- 
ral thresholds ought to be found in other sense modalities. 


The auditory neural unit, then, is assumed to have a threshold 
intensity Ip which varies at random in normal fashion. If, during 
the instant of stimulation, the stimulus intensity is at or above the 
threshold, (i.e., for all J > Ir) the neural unit fires. If it is below 
the threshold, the unit does not fire. Our assumptions indicate that 
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the probability of response to an intensity [> is the probability that 
I, is above the threshold of the neural unit and this is: 


log Ig ] [— 2 
l : (=) eGordy: (9) 
o V2n o 


-0O 


pily = 


Equations (8) and (9) give us RA as a function of J and of para- 
meters m and o. This function is represented by the curve of Fig- 
ure 1 with m = —25 db, « = 20 db. For comparison the experimen- 
tal points are indicated by dots. It should be noted that p(J) is a 
function of log intensity in our definition. This formulation is im- 
posed by our data. On the other hand the curves of Pecher (1938) 
and E. T. von Briicke et al. (1941) appear to be on linear intensity 
scales. We are prepared to offer no simple explanation of the dis- 
crepancy, other than to point out that our measures of intensity 
refer to the acoustic click and not to proximal electrical stimulation. 
A number of alternative hypotheses may be offered to account for 
the log transformation. We see no reason for committing ourselves 
to any one at the present time. 


39) 


90 
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Figure 2, The reconstructed intensity growth function of N,. The ordi- 
nate is a probability grid on which a normal integral is linear. The experimental 
points are taken from Figure 1. Similar functions have been obtained for dif- 
ferent intensities of the second click from this same animal, as well as from other 
experimental animals. 


Figure 2 shows p(J) reconstructed from the curve of relative 
amplitude in Figure 1. The points are plotted on a probability grid 
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on which the normal integral is rectilinear. In all it appears that 
our assumptions are reasonably validated by our results. We have 
fitted the data of Figure 2 to the intensity growth curve of the first 
neural component. The results are as expected (Rosenblith and Mc- 
Gill, unpublished data); N, follows the theoretical curve until the 
microphonic begins to grow, after which the neural component oscil- 
lates around the theoretical value as intensity increases. 


Discussion 


We have assumed that each neural unit fires with nearly the 
same average voltage. It is quite unlikely, however, that the magni- 
tude of the response emitted by any single neural unit is constant or 
that any two units contribute exactly the same amount to the total 
response. We can, on the other hand, interpret the response magni- 
tude of a neural unit in the same way as we have dealt with its excit- 
ability, i.e., by stating the probability of a response of magnitude 
Ry, namely p(R,.). The mean value of response magnitude will then 
be, assuming a finite expectation: 


te = E(R;). (10) 


If the probability of firing is statistically independent of the prob- 
ability of a given response magnitude, then the expected magnitude 
of the total response in a population of N neural units is N-p(Z) - ue. 
Independence here implies an all-or-none law for neural activity, 
which is somewhat unusual in comparison with the common formu- 
lation, but it is nevertheless effective. Actually an all-or-none law of 
sorts has been implicit in our entire development. The refractory 
period operator p*(t) has been described as a function of two vari- 
ables: time separation of the two clicks and intensity of the second 
click.* In the two-click problem we concern ourselves not with the 
stimulus intensity that fires a single neural unit, but only with the 
fact that it fires. We cannot, however, neglect the intensity of the 
second stimulus. Our experiments show clearly that a given stage of 
recovery of the second click response is reached much more rapidly 
as the intensity of the second click increases. 

The fact that component N, of the click response does not be- 
come supernormal when paired with a preceding click permits us to 
treat p*(t) as a simple recovery probability, or geometrically as 


a dis- 


*The term refractory period is employed here in a generic sense. ae H 


cussion of refractoriness, depression, fatigue, and their interrelations, 
Lorente de N6 (1947, p. 367). 
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a probability plane with time and intensity of the second click as 
coordinates. Undoubtedly the most general formulation of p* (¢) is 
in terms of a transform on the parameters (m and c) of the second 
click response. Nevertheless, we have chosen to work with the re- 
covery probability in order to preserve the simplicity inherent in 
responses to pairs of clicks. In the present paper we have not been 
concerned with the mathematical expression for the p*(t) operator 
(which has been absorbed into the constant ¢). The operator is, how- 
ever, amenable to direct experimental analysis and will be dealt with 
empirically in a subsequent paper. 


We do not require that each neural unit have the same mean 
threshold intensity or the same sigma. As a matter of fact we can 
conceive K classes of neural units with different threshold distribu- 
tions. Each class might contain an appropriate number, ,, of neu- 
ral units having an average response magnitude, uz,. In this case 
the number of elements and the average response magnitude consti- 
tute weighting factors for each threshold distribution. We rely upon 
a well-known statistical theorem which states that a linear combina- 
tion of normal variables is itself normally distributed, to guarantee 
that the average of classes of neural units behaves as our experi- 
ments seem to indicate. 


So far we have dealt with a single problem: the amplitude of 
the first neural component of the round window response to clicks. 
Thus, we have considered only one aspect of a complex response, dis- 
regarding in our present formulation changes in latency, temporal 
dispersion, and behavior of other neural components. We feel rather 
confident that the present model does not contain features that would 
prevent us from extending our predictions to these other dimensions 
of the response. It is, as a matter of fact, not too hard to see how the 
flexibility of the p*(t) operator might permit us to investigate such 
phenomena as fatigue, masking, or recruitment, or to account for 
the effect of certain physical variables such as temperature upon the 
click response. There are many ways in which we can manipulate 
the experimental situation to get at the p(t) operator and its para- 
meters—and this, after all, is at least one of the functions of a mathe- 
matical model in electrophysiological research. 
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A SUGGESTION FOR A MATHEMATICAL APPROACH TO THE 
PROBLEM OF THE EFFECT OF CLIMATIC 
FACTORS ON HEALTH 


N. RASHEVSKY 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


It is suggested that the effect of sudden temperature changes upon 
the organism may be due to the variation of two factors, whose rates 
of changes are functions of the temperature. It is assumed that the 
variations of the two factors follow differential equations which are for- 
mally identical with those of the two-factor theory of excitation, and 
in which the temperature formally plays the role of the stimulus. By 
proper choice of the parameters the two factors are made equal for anv 
constant temperature. Assuming that the greater the absolute value of 
the difference of the two factors, the greater is the instability of the 
organism, and the greater is its susceptibility to disease, relations be- 
tween the suddenness of temperature changes and disease incidence are 
derived. The incidence of weather-affected disease for a given period of 
time is found to be a linear function of the mean hourly temperature 
variation during that period, a relation which can be in principle veri- 
fied. 

Paraphrasing the well known remark of Mark Twain we may 
say that a great deal has been written about effects of climate and 
weather on health, but very little, if any, scientific study has been 
undertaken. The most prolific writer on the subject in this country 
was undoubtedly the late William F. Petersen. The several huge vol- 
umes of his The Patient and The Weather (Petersen, 1936-38) con- 
tain an enormous number of diagrams and charts, which, unfortu- 
nately, could not convince a scientist. While they purport to show a 
close parallelism between climatic variations and variations of health, 
not a single correlation coefficient is computed, and no quantitative 
treatment of the data is given. In his Man, Weather, Sun (1947) 
Petersen does give correlation coefficients between climatic changes 
and some pathological changes, but a recomputation of some of the 
coefficients, taken at random from the data given in text, gives quite 
different values than those claimed by Petersen. 

A much more modest approach to the problem was made by 


Manuel Morales and Eugenia Tarver (1948), who observed the in- 
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cidence of some diseases among the crew of the battleship Washing- 
ton as the ship traveled through different climatic environments. 
The results show a very clear increase in incidence of respiratory 
disease associated with sudden drops in temperature and vapor pres- 
sure. There was no attempt, however, to establish a quantitative 
relation between the incidence of disease and the climatic variations. 

One of the reasons for such a lack of more definite results lies, 
undoubtedly, in the lack of any adequate theoretical background for 
planning the necessary observations or experiments. It is the pur- 
pose of this note to suggest a working hypothesis which, when de- 
veloped into a theory, will suggest some quantitatively verifiable re- 
lations. 

Available qualitative observations indicate that the incidence of 
some diseases is related to sudden changes of climatic variables such 
as temperature, regardless of the direction of the change. For def- 
initeness we shall discuss here only the temperature as a climatic 
variable. Mutatis mutandis, this reasoning may be applied to other 
variables. 

A rather superficial analogy to the phenomena of nerve excita- 
tion under the action of sudden changes in electric currents suggests 
the possibility of a formal application of a “two-factor theory”’ to 
our problem. 

Suppose that there are two physiological factors, a and b , which 
determine the normal functioning of some parts of our organism. Let 
the organism function normally when the two factors are in a given 
constant ratio, which, without loss of generality, may be assumed to 
be unity. Hence, for normal functioning, we require 


a=b. (1) 


Whenever a # b, the organism is not functioning properly, its 
resistance to disease is lowered, and the incidence of disease increases 
with increasing absolute value ja — b|. At a first approximation 
we may consider the incidence Idt of a given disease during the in- 


terval dt of time to be proportional to |a — b\dt, so that, with a as 
a coefficient, 


Idt=a\a— bldt. : (2) 
Denoting by T the temperature and by K , M,k, and m four 


constants, let the variations of a and b with time Ba given by the 
following differential equations: 
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ae RE eT b 
a =k, —= = 
dt dt ie (3) 
with 
K M 
—_—_— — — ; 4. 
iia es (4) 


If the temperature is kept at a constant value T for a sufficiently 
long time, a and b will reach their asymptotic values: 


K 
a=b=—T=—T=xT. (5) 
k m 
If, at the moment ¢ = 0, the temperature is suddenly changed 
by an amount T, 2 0, now becoming equal to 7 + 7,, then, from 
that moment on, a and b will vary according to the equations 


a=«l + «T,(1—e*) (6) 
and 
b=cT i+ «7, (1l1—e”"). (7) 
Hence 
|a — b| =x|T,(e™ — e**) |. (8) 


As t tends to infinity |a — b| tends to zero. As is readily seen 
from (8) the quantity |a — b| has a maximum for 


1 m 
f= log —> 0. (9) 
m—k kL 


Now let a continuous change in temperature be represented by 
Del) 1 Ge 67), (10) 


so that the total temperature now is 7 + T*(t) where qT, and 4 are 
constants and 7,* 2 0. Introducing (10) into (8) and integrating 
with the initial condition 

a=b=xT fort=0, (11) 
we find 


* 


KT, 
Genie tel i= ex) aah (61 —¢*") (12) 


and 


* 


MT, 
Dokl Kia. (Le )o F (e4*—e™*). (13) 
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If 
Fee for (14) 
then | 
Ey ey at (15) 
ca 


Hence in this case the last term in equation (11) is small com- 
pared to the preceding one. As 4 increases that term decreases. The 
same holds true of the last term of (12). For sufficiently large val- 
ues of 2, equations (12) and (13) reduce correspondingly to (6) and 
(7), as should be the case physically. For sufficiently slow varia- 
tions of 7, when 4 is not too large, equations (12) and (13) may be 
conveniently used. 

If, on the other hand, 4 is very small, then (12) and (18) reduce 
correspondingly to: 

Oe +eTy A—e"); (16) 


b=«l +xT (1 —e); (17) 


so that a= b all the time. A sufficiently slow variation of T does not 
result in a physiological disturbance nor in an increased susceptibil- 
ity to disease. 

From equation (2) it follows that total incidence of disease dur- 
ing the time following a sudden change in temperature will be pro- 
portional to 


fe ja — b| dt. (18) 


It follows from (8) that the integral (18) in its turn is proportional 
to «xT, or to «A, if we now denote the sudden change in temperature 
by A. For practical purposes the upper limit of integration may be 
taken as 1/nl, where | is the smaller of the two quantities, k and m, 
and n is of the order of 3 or 4. 

Absolutely sudden changes do not actually occur. But, if we 
make a plausible guess that a ~ b ~ day-?, then we may consider a 
strong hourly variation as sufficiently sudden. Let N(A)dA be the 
distribution function of the absolute values of hourly A’s which have 
been observed over a sufficiently long period, so that N(A)dA denotes 
the number of variations whose value is between A and 4 + dA, re- 
gardless of the sign of 4. We may now write the following for the 
total incidence J of disease during the period: 


f=0 { AN (A)dA~=aA, (19) 
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where A denotes the average hourly variation during the period. 

Hence, under the above assumption, as expressed by equations 
(2), (3), and (4), we find a direct proportionality between the in- 
cidence of disease and the average hourly variation 4. Both can be 
measured directly. 

The procedure thus indicated by the theory is to look not for 
individual temperature variations and correlate them with individ- 
ual cases of incidence of disease, but to determine the incidence of 
disease for various periods of the year, for which periods the aver- 
age hourly variations of temperature are sufficiently different, and 
to see whether the predicted linear relation holds. 

If we assume that the incidence Jdt is not proportional to 
|a — b| dt but is some more general function 


Idt = F (a — b) dt (20) 


then, by introducing (6) and (7) into (18) and integrating with 
respect to ¢ from 0 to o , we find the incidence per sudden change 
as a function U(A). Then, instead of (19), we find 


I=aJU(A). (21) 
The function U(A) can be determined from the empirically ob- 
served N (A). 
We may make a more general assumption, namely, 
Idt=a[|a — b| —h] dt, (22) 


where h is a threshold, so that Jdt = 0 for |a — b| < h. Such a rela- 
tion is particularly likely to hold for the incidence of death, though 
it may hold also for the incidence of disease. The total incidence 
I(A) which results from a sudden change A is now given by 


1(4)= {°[la—b|—A] at, ac) 


where ¢, and ¢, are the roots of the equation 
la(t) —b(t)|=h. (24) 


From (6) and (7) it follows that equation (24) has roots only 
when A > h/n, where 7 is the maximum value of «le — e*|, ob- 
tained by introducing (9) into it. Hence [(A) is also a function 
of h, so that 
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l[=ViA, Rh): (25) 


Consequently, we have for a long period 
lV hk): (26) 
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A distance bias is imposed on the probability of direct connection 
between every pair of points in a random net. The probability that 
there exists a path from a given point in the net to another point is 
now a function of both the axone density and the distance between the 
points. A recursion formula is derived in terms of which this prob- 
ability can be computed. 

The rate of spread of an epidemic where probability of contact de- 
pends on the distance between the individuals can also be computed from 
the recursion formula. 


Using an approximation method, R. Solomonoff and A. Rapoport 
(1951, see pp. 107-17) have deduced an equation which implicitly de- 
fines y , the weak connectivity of a random net as a function of its 
axone density, i.e., the average number of axones issuing from each 
point. It was shown how the weak connectivity was a measure of the 
probability that an arbitrary individual in a closed population suc- 
cumbs to a contagious disease introduced by a single infected indi- 
vidual, where it is assumed that the disease has a finite contagious 
period, followed by immunity or death, and where contacts between 
any pair of individuals are assumed equiprobable. 

To make the epidemic picture more realistic, it is obviously nec- 
essary to introduce some sort of “distance bias,” i.e., a distribution 
governing the probability that a given individual will come in con- 
tact with another given individual in the population. Such prob- 
ability cannot in general be assumed constant, since individuals are 
more likely to come in contact with some individuals than with 
others. A constant probability of encounter represents the extreme 
case of very thorough “mixing.” 

Assuming the approximation method of Solomonoff and Rapo- 
port valid, we shall derive a procedure for computing 7 (r) , the prob- 
ability that an individual at a distance r from the original source of 
infection will eventually contract the disease, if the probability of 
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contact of two individuals is governed by a function of the distance 
between them. 

In terms of a neural net, y(7), will indicate the probability that 
there exists a path from a neuron to another neuron at a distance 7. 

In order to generalize the result of the above-mentioned paper 
(Solomonoff and Rapoport, 1951), we will give another derivation of 
equation (20) of that paper, which, incidentally, will enable us to 
compute the rate of spread of the epidemic (or, in terms of the neural 
net, the expected number of neurons ¢ axones removed from the giv- 
en one). 

We will again follow the axone-tracing procedure described in 
the above paper. Let p(t) designate the probability that a given 
neuron is contacted on the tth tracing (whether or not it had been 
contacted on previous tracings). Note that in the absence of distance 
bias, the »’s represent independent events, since being contacted on 
one tracing does not affect the chances of being contacted on another 
tracing. Likewise if g(t) = 1— p(t), the q’s are independent prob- 
abilities. 

Then the probability that a neuron is contacted for the first time 
on the tth tracing will be 


t-1 t 


p(t) Hai) =U—a()] Hai). (1) 


Since each neuron sends on the average a axones, and there are N 
neurons in the net, the expected number of axones to be traced on 
the (¢ + 1)th tracing will be 


aN[1—a(¢)] Ia (i). (2) 


The probability that any neuron in the aggregate is not contacted: 
by any of these axones on the (¢ + 1)th tracing will then be 


t-1 
aN [1-9(t)] II a(4) 
i=0 


q(é+1)=(1—1/N) ' (3) 
which for large N can be written as 
a(t +1) = Exp {a 1— a aw } 
t-1 t 
= Exp { ~a | Za —TaW ts (4) 


be Taking a product of both sides of (4) with respect to t, we ob- 
ain 
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t+1 t 


a3) = I Bx | —a [ Haw —Matiy | 


=Exp | Sy ia) ~1a(i) | (5) 


= Exp | —a | 900) — 140i) |t- 


But q(0) = (N —1)/N ~ 1 for large N.. Therefore, since 


t+1 t 


Lim I q(j) =Lim II q(j)=1—», (6) 
t>o j=l t>o j=0 


we have, taking the limit of both sides of (5), as t approaches in- 
finity, 
l1—y=e” or y=1l—e™”, (7) 


which is equation (20) of the above-mentioned paper. 

Let us now introduce a probability density f(7), such that 
f (7) dS will denote the probability that an axone will contact a neu- 
ron in the vicinity of a point at a distance r from the neuron which 
sent out the axone. Here dS is a differential of length, area, or vol- 
ume, depending on the dimensionality of our net. Evidently we must 
have 


if FPS, (8) 


where the integration is taken over the space occupied by the net. 

To fix ideas, let our net be one-dimensional and let us choose the 
origin at the point where the tracing of neurons is to start. Then 
f\|x — x,|Axz will denote the probability that an axone issuing from 
the neuron at x, will synapse on a neuron in the vicinity of the point 
x, or vice versa. 

We wish to extend the method which we used to derive y(a) as 
a function of axone density a, to a derivation of y(a,”), or, since a 
is fixed, of (x). That is to say, we seek the probability that there 
exists a path through any number of internuncials from a neuron 
taken at the origin to a neuron in the vicinity of «, where distance 
bias has been imposed by the function f|”z — | . 

However, the method used above is not immediately generaliz- 
able to the case with distance bias, because when distance bias is 
introduced, the p’s and hence the q’s are no longer independent. To 
show this, assume there is a bias in favor of the nearer neurons 
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being contacted on each tracing (i.e., bias in favor of shorter axones). 
Then if the probability of being contacted on a particular tracing is 
high, this means that many neurons in the immediate vicinity will 
be contacted, some of them for the first time. And this in turn means 
that many axones will be traced from the immediate vicinity on the 
next tracing. Hence being contacted on some tracing enhances the 
probability of being contacted on subsequent tracings. 

To avoid this difficulty, we shall rewrite equation (3) in terms 
of another probability P(t), which will be defined as the probability 
of being contacted for the first time on the tth tracing. By equation 
(1) we have 


t-1 
P(t) = Mra ia). (9) 


It also follows from the definition of P(t) that 
1-3 PC) | a(t +1) =1—S PU), (10) 
aét+y=|1-EPQ) | 1-sPewH |, aw 


t -1 
1—q(t¢+1)—=P(é+1) | Eee) | , (12) 
j=0 
Rewriting equation (4) in terms of P’s, we then obtain 
t 
P(t+1) =| 1—3 P(j) | (1— 6-0, (13) 
4=0 
Equation (13) is a recursion formula for the P’s which allows 


the calculation of all successive P’s once the first is known. The total 
number of neurons contacted by the tth tracing is evidently 


t 
NDP). 
g7=0 
In terms of the epidemic this function measures the spread of the 


epidemic (total number who have succumbed at time t). Further- 
more we have 


ZPO)=y. (14) 


In contrast to equation (4), equation (13) is generalizable to 
the case with distance bias. Let p be the linear density of neurons 
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in our net. Subdivide the net into “macroscopically” small regions 
Ax, analogous to those considered in the kinetic theory of gases. 
That is, the regions are small compared to the extent of the net but 
each region contains a great number of neurons. Then the expected 
number of neurons contacted on the tth tracing for the first time 
in the region Ax about the point x will be pP (x,t) Az. Consequently 
the expected number of axones to be traced from that region will be 
apP (x,t) Ax. The probability that a neuron at x is not contacted 
by any axone issuing from any region on the (¢ + 1)th tracing will 
be 


I 1—f\x—a|Aa (15) 


ee 
where the product is taken over all the regions Ax , taking a repre- 
sentative point x from each region. Therefore the probability that 
a neuron at x, will be contacted by at least one axone on the (t + 1)th 
tracing will be 


ussitl Pe fa ae (16) 


apP(a,t) Ag 
& 


Now the probability that the neuron at 2, had not been con- 
tacted on a previous tracing is 


1-3 PW). (17) 


Furthermore the events represented by the probabilities (16) and 
(17) are now independent, since the distribution of axones to be 
traced is now given in terms of the distance bias. Therefore the 
product of the probabilities (16) and (17) gives us 


PoC, t+. 1) 
t apP(2,t)Ag 
= | ZW) |[a-a[ t-te alae | |, (18) 
j=0 @ 


which is the generalization of (13) desired. In the case without dis- 
tance bias, where pAxz = 1 and f|x — %|4x = 1/N , equation (18) 
reduces to (18), as it should. 

Equation (18) is a recurrence formula for the P(x, t). Hav- 
ing calculated the successive P(x, t), we obtain y(x) from the re- 
lation 


(0) =B Pa, t). (19) 
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Critique of the Method. The above method fails if the density of 
neurons is too low, and/or the distance bias is too strong. These two 
parameters have a definite relation to each other. If, for example, 
the distance bias is weak, that is the probability of direct connection 
is fairly constant over considerable intervals, then we may artifi- 
cially increase the density of the neurons by taking Aw larger. But if 
both the density is low and the bias strong, we cannot take larger 
intervals without committing gross error, because in our treatment 
the probability of direct contact was considered constant through- 
out an interval. 

To illustrate this we shall take an extreme case, namely where 
direct contact is possible only between two adjacent neurons. Then 
we are forced to take Az so that it contains only one neuron. Start- 
ing our tracing procedure from the origin, we see that successive 
neurons to the right of the origin will be contacted on successive trac- 
ings only as long as not all axones of a neuron will go to the left. The 
probability that this event will not occur is evidently (1 — 2), and 
the probability that it will not occur ¢ successive times is (1 — 27)*. 
Hence (1 — 2“) ig the probability that a path exists to exactly t 
neurons to the right of the origin. The expected number of neurons 
to which a path exists will then be given by 


N 


Peer 20 Bi eA (20) 


t=0 


counting both the neurons to the right and to the left of the origin. 
For N very large, the expected number of neurons to which a path 
exists remains finite, namely, 


poy res ee | 
's bee 21 
2% 420 2° on (eA) = 


where A = (1 — 2). This gives for the expected number of neurons 


2 Le 
Lt \a2@ 420-1). 


28 2-21 


It can be shown that this result cannot be obtained from equa- 
tion (18) by setting p= Az = 1 and putting for f |% — ao| the value 
1—2~ for « — a =1 and zero otherwise. The reason for this lies 


*For small values of a it is important to correct the quantity on the right 
fi of (22) by subtracting unity, since the neuron at the delet was counted 
ce. 
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in the fact that in our treatment it was supposed that at each trac- 
ing a number of neurons were contacted in each region and that 
probability of contact was constant throughout a small region. In 
the extreme case just treated, however, a neuron was either contacted 
or not, and the neuron was coextensive with the region. If a neuron 
was contacted, the same number of axones were traced on the next 
step; if it was not contacted, the tracing stopped altogether. The dis- 
crepancy is analogous to the one observed in the theory of radiation 
of a black body, which made necessary the postulate of discontinuous 
radiation. So long as the density of the neurons is sufficiently great 
and/or the distance bias sufficiently weak, the connectivity problem 
can be treated as a continuous case. 

Our model of the net with distance bias actually assumes a sub- 
division of the net (or the community in the contagion model) into a 
number of subnets (or subcommunities) such that within each sub- 
community contacts are equiprobable, but intercommunity contact is 
governed by a distance bias. The picture is not too far-fetched if the 
subcommunities are regarded as “households” where thorough mix- 
ing of contact occurs. 

This investigation is part of the work done under Contract No. 
AF 19(122)-161 between the U. S. Air Force Cambridge Research 
Laboratories and The University of Chicago. 
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Two theories are discussed to account for the observed reso- 
nances of biological cells at sonic frequencies. One theory assumes the 
cell wall to be a stretched balloon surrounded by, and filled with, an in- 
compressible fluid. The other treats the cell wall as a rigid shell. Both 
lead to reasonable physical constants for the cell wall. 


Introduction. The biological effects of intense sound waves in 
liquids have been a subject of interest for many years. Aqueous sus- 
pensions are dramatically altered by brief exposures to intense sonic 
fields; protozoans, algae, and bacteria are destroyed; red blood cor- 
puscles are torn into fragments; and small animals are killed (Cham- 
bers and Harvey, 1931). In almost all experiments the occurrence 
of cavities in the liquid has accompanied the destructive effects. There 
seems little doubt that the destructive results of intense acoustic fields 
in liquids are produced by the shearing forces in the immediate 
neighborhood of these cavities. 


Qualitatively, the same results are obtained at 1 kc./sec. and at 
1 mc./sec. (Chambers and Gaines, 1932). Most studies have been 
carried out at an isolated set of frequencies, and little attempt has 
been made to measure the intensities of the acoustic field. A study, 
reported in another journal, shows that there are optimum frequen- 
cies at which the breakdown occurs much more rapidly for a fixed in- 
tensity. These frequencies depend on the type and size of the or- 


*The initial phases of this work were supported by a grant from the Wis- 
consin Alumni Research Foundation to Professor H. B. Wahlin and were in- 
cluded in the author’s theses submitted to the Graduate School of the Univer- 
sity of Wisconsin in partial fulfillment of the degree of Doctor of Philosophy, 
August, 1949. The work is being supported at present by a grant to the Johnson 
Foundation from the Raytheon Manufacturing Company. 
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ganism treated. Table I shows the physical sizes and the optimum 
breakdown frequencies for several strains of paramecium. 


TABLE I 
SPECIES 2a 26 a eff. a/h » obs. 
P. caudatum 222.6 u 63.04 71.4p 8.3 1.2 ke./sec. 
P. busaria A720 B14 41.6 9.8 1.7 
P. calkensi 125.6 56.2 45.4 6.2 1.9 
P. aurelia G’s 123.8 Zoe 39.8 6.9 3.0 
P. aurelia 51 127.8 30.2 38.2 6.8 3.5 
P. aurelia 8la 98 38 34.0 5.0 4a 
P. trichium 79.6 Sys 29.2 — Toes 


Physical sizes and resonant frequencies of seven strains of Paramecium 
species. 2a is the longest diameter; 2b is the longest diameter at right angles 
to 2a. The effective radius a equals a+0/2; h is the thickness of the cell cor- 
tex. All figures are the average sizes of at least twenty individuals. 


The wave length of sound in water at 1 kc./sec. is 1.5 meters, 
i.e, about 10* times larger than the average P. caudatum cell. Thus 
the resonances producing increased breakdown must be due to prop- 
erties of the cell itself. As the breakdown does not occur in the ab- 
sence of cavitation, and as the cavities occur most readily outside the 
cell (Harvey, 1947), it seems reasonable to assign these resonances 
to the properties of the cell membrane. Hence these resonances offer 
a method of investigating the elastic properties of the cell membrane 
(or possibly the cell cortex). 

To actually assign any significance to these optimum breakdown 
frequencies, the cell must be interpreted in terms of a simpler physi- 
cal model. The present paper deals with the mathematical develop- 
ment of two such models, both of which make resonances reasonable 
in the observed frequency range. The first model treats the cell 
membrane as one possessing an interfacial tension, but no rigidity. 
This model has been used in many previous studies (Harvey, 1938), 
although the methods used could not be applied to the cells consid- 
ered here. The second cell model to be considered assigns a rigidity 
to the cell cortex. In both models the cells are considered to be spheri- 
cally symmetric, to be filled with an ideal, incompressible liquid, and 
to be surrounded by another similar liquid. Clearly the paramecia 
do not have spherical symmetry, and the viscosity of the cell con- 


tents is certainly important. Thus the present development is only 
a first approximation. 
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Interfacial Tension Waves. As noted above, this first model re- 
places the cell by a spherical shell, lacking any rigidity, but possess- 
ing an interfacial tension. It makes no difference if this tension is 
a true liquid-liquid interfacial tension (as that between water and 
oil), or a liquid-membrane tension, or a tension residual in a 
stretched membrane (such as a rubber balloon). Physically all of 
these may exist at the cell boundary. Values of this interfacial ten- 
sion, T,* measured on large non-mobile cells, range from 0.01 to 3.0 
dynes/cm. The theory discussed here gives values of 7 from 3 to 
10 dynes/cm. for paramecia. 


The shapes of the paramecia are far too complex to handle 
mathematically. Even the simpler model of a prolate spheroid pre- 
sents real difficulties. Rather than become lost in the details of a 
closer approximation, it seems wise to consider the simple case of a 
spherical cell. In the case of the electromagnetic wave equations, 
this symmetry approximation does not alter the order of magnitude 
of the sizes of resonant conductors. It appears reasonable to assume 
that the same thing would be true here. 


As the compressibility of the liquids is ignored in this treat- 
ment, the starting point is not the acoustic wave equations which 
depend on the compressibility of the liquids for their existence. In 
any event these acoustic waves are many orders of magnitude too 
large. Surface tension waves of a liquid have a much shorter wave 
length than acoustic waves. The discussion of interfacial tension 
waves is similar to that of surface tension waves. The appropriate 
equations of motion for an ideal incompressible liquid can be found 
in many texts. The notation used here is similar to L. Page’s (1935). 

The liquid motions are assumed to be irrational and only rela- 
tively small velocities are permitted. Page shows that these motions 
can be described by the following equations: 


v= — V9; (1) 
V76=0; (2) 
b=p/p+Q+f l(t). (3) 


In these v is the vector velocity, & the velocity potential, p the pres- 
sure, p the density, 2 the potential energy per unit mass, and f (t) is 
an arbitrary function of ¢, the time. Since only the spatial deriva- 
tives of ® have any physical significance, we choose 


*For a complete list of symbols, see pages 104-05. 
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f(t) =0. 

The only potential energy in this problem is that due to gravity, i.e., 
Q=gh. 


Since this is approximately constant over a small cell, it may be ig- 
nored in a problem involving alternating pressures and velocities. 


Quantities inside the cell will be designated by the subscript 1, 
those outside by the subscript 0. Standard spherical coordinates r, 
6, and y will be used. All equations, from this point on, will refer 
to the time dependent portions of v, ®, and p. In particular, equation 
(3), the equation of motion, can be rewritten as 


®; — Di/ pi 
(4) 
by = Do/po 


The above are general hydrodynamic relationships. Now consider 
the specific cell model. The liquids, being ideal, may slip freely over 
the membrane but not lose contact with it. Let a be the radius of 
the equilibrium position of the spherical membrane, possessing the 
interfacial tension 7. Assume a small deformation which carries 
the point (2,6, w) to (a+ R, 6, wy), where R, the radial extension, 
depends on 6 and w. The excess force per unit area on the membrane 
in the +, direction is 


Jk O tine 0 de Lae goth 
r= ace men | EST SH 1h Ge 2 
a’? sin 6 00 ( =) a? sin? 0 0 w’ 


(5) 


Hence we find 
Di SS Po =F F ° 


Other boundary conditions are: a) the center of the cell stands still; 
b) no effects of the motion are observed as r > 3; ¢) the liquid 
adheres to the membrane at all points and at all times. This last 
condition will be satisfied if the liquid radial velocity is continuous 
across the membrane. Expressed analytically these conditions are: 


EUGENE ACKERMAN 97 


a) ea=40 at p= 0; (6) 
or 

b) d,—> 0 as ete (7) 
0 ®; 0 Bo 

®) perks edligyy oY Wiqiees (8) 


A displacement of the membrane in the 6 or y direction would 
contribute nothing since the membrane lacks rigidity. As an added 
simplification for this model assume that #;, &, and R are inde- 
pendent of y. 


An appropriate solution for equations (2), (6), and (7) is: 


&;=> A, 7” e/*' P,, (cos 6), 
n=0 
(9) 
=> Brar'™» eft P, (cos 6), 
n=0 
where A,, B,, and w, are constants and P,(cos @) is the nth Le- 


gendre polynomial of argument cos @. The constant w, is 2a times 
the characteristic frequency of the nth mode. 


Substituting equation (9) into (8) gives us 


B,=— Go As: (10) 


Oa ih 


We now combine (4) and (5) to eliminate the pressures and get 


: Do st 0 : go 
oh) fh 1 Sin 
bases pi a pi Sind 00 06 
Po * tM 0 k oR 
—— >, — —_——_ —__|{_ sin 0 -—— ° 
pi a p; sin 600 006 


Differentiating this with respect to time, and using (1) to eliminate 
R, we get 


Po ” T 0 ‘ 0? ®; 
6, = — 6, — ———_ —_[ sin@ : 
pi a? p; Sin # 08 aroé 


To eliminate the velocity potentials, we substitute (9) and then (10) 
into the last equation as follows: 
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— w,? A, a" =— o,? Bn a — nr? (n + 1) a" An, 


3 


pi 


n Po fh n 3 1 
==) (en Soe | eens ee or I ae Ae 
( (is an Oy: pi a nN 


In a resonant vibration, one frequency, w,, will dominate the mo- 
tion; only one value of n will be important. The values 0 and 1 for 
n are prohibited by equation (6). The lowest resonant frequency 
corresponds to n = 2; it is the solution of: 


or ( a i 3/2 )(2 ie ). (11) 
pi a 3 pi 


To apply this to paramecia, we choose pi = po = 1 gm./cc. and find 


Wy.” = Te T ae 


The —3/2 power relationship seems unusual, but for surface tension 
waves on water of frequencies near 1 kc./sec. we have 


eT nea 
o? = a t 
p ( A 
where 4 is the wave length and T, is the surface tension. (For any 
given wave velocity, there are two wave lengths possible. Only the 


shorter one corresponds to real waves at frequencies of the order of 
magnitude considered here. ) 


The above formula for w, may be used to calculate the values of 
T for the various paramecium strains. The results assign a value of 
T = 3 dynes/cm. to P. caudatum, and T = 10 dynes/cm. to P. 
trichium. It was originally hoped that these would be closer to one 
another, but there is no a priori reason to assume that all species, or 


even all strains of the same species, should have the same interfacial 
tension, 7’. 


Rigid Shell Shear Waves. The experiments measuring inter- 
facial tensions on non-mobile, or slowly moving, cells can be inter- 
preted in ways other than the usual one; some as measuring the 
tensile strength of the membrane and others as measuring the rigid- 
ity of the cell cortex. The optimum breakdown frequencies of para- 
mecia can also be interpreted as resonant vibrations of a rigid spheri- 
cal shell, emersed in, and filled with, incompressible liquids. Several 
modes of vibration will be considered for this model. 
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Lord Rayleigh (1890), points out in his chapter on rigid shells, 
that there are, in general, two types of vibrations possible. The more 
common type involves no extension of the mid-surface; it is the mode 
usually discussed for flat plates and bells. However, according to a 
theorem due to Jellet, no closed convex surface can be deformed in 
any way without extension (or compression) of the mid-surface. 
Therefore, the usual methods of treating shells are valueless here. 

For vibrations with extension of the mid-surface, both the ki- 
netic and the potential energies are proportional to the shell thick- 
ness, h. Closed (and almost closed) shells vibrating in air will, there- 
fore, have resonant frequencies which are independent of h. Ina 
cell such as a paramecium it is natural to interpret h as the thick- 
ness of the cortical layer. 

The rigidity of the cell cortex is negligible compared to steel, 
glass, or even wood. However, protein gels do have a finite rigidity. 
Values of the coefficient of rigidity, 1, have been measured for fibrin 
gels (Ferry and Morrison, 1947). The modes to be discussed here 
lead to values of uw in the same range as that for fibrin gels, i.e., 10° 
to 10° dynes/cm?. 

The analysis follows that of H. Lamb (1882), who fully dis- 
cussed the vibrations of spherical shells in air. Only modes symmet- 
rical about a diameter will be considered. These fall into two types: 
1) those involving tangential motion only; and 2) the modes with 
both radial and tangential motion. 

1) The tangential type of modes will not be affected at all by 
the intra- and extra-cellular liquids, since these are assumed to be 
ideal, and hence to slip freely over the cell wall. In this type of mo- 
tion a point (a, 6, y) is carried to (a,6,y+ ¥). The angular dis- 
placement in the yw direction is ¥; it is time dependent. Since the 
liquids play no role, Lamb’s analysis for shells in air applies directly 
to the cell model. He shows that 


d 
w= A, ———— [Pn(cos @)] ef" 
d (cos @) 


and 
Nad 
Wy? = —_ (n—1) (n +2). 
ps a” 

The subscript s refers to the shell, i.e., p, is the density of the spell 
This Y mode gives for P. caudatum for = 2, u oe 0.8 X 10 
dynes/cm?. This is in order of magnitude agreement with the fibrin 
gels; a possible mode. 
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2) The modes with radial motion as well as tangential involve 
both the motion of the liquids and that of the shell. A more detailed 
description of these modes is therefore necessary; it occupies most 
of the rest of the paper. 

Many of the equations from the section on interfacial tension 
waves may be used without change. The useful ones are (1), (2), 
(3), (4), (6), (7), and (9). These apply to the liquid motion, and 
additional equations are necessary to describe the motion of the shell. 
Changes occur in 7 and 6; these are denoted by RF and ©. Strains 
will occur in the r, 6, and y directions, and are denoted by ¢, é, 
and es; respectively. The average stress on the shell in the +r direc- 
tion is: 


— p,=—1/2(pi + Po). (12) 
The accelerating force per unit shell area in the +,r direction is: 
Ap, = Pi —Do- (13) 


These last two equations describe analytically the difference between 
the cell model which is filled with, and surrounded by, incompressible 
fluids and a shell vibrating in air. In the latter case both —p, and 
Ap, must vanish. 

The equation used to describe the adherence of the liquids to the 
membrane in the interfacial tension model, (8), must be modified to: 


0 9; 0 By 
or or 


Sete 


It will be assumed that |he,| << |R,|; therefore, equations (8) and 
(10) may be used for this model also. This assumption will be con- 
sidered at the close of the development; it will rule out one of the 
solutions as self-inconsistent. 


From the geometry of the system (Lamb, 1882), it follows that 
the strains can be represented by: 


R./a+ Gay 
é. = R,/a + — 14 
ra (14) 
and ! 
é;—R,/a+ 0 coté. (15) 
The above equations describe the motion of the liquids, which 


are similar to those in the interfacial tension model, and also the 
stresses and strains on the rigid shell. To proceed further the equa- 
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tions of motion of the shell are needed. These can be derived from 
Hamilton’s principle as follows: 


at+h/2 a+h/2 
< OEP, 
hph,8 KR, 42°dr= here em AR 8R,4ardr 
a-h/2 a-h/2 ORs ; ; 


therefore 


8 


* Ms) 
| asian el Ce + A De; (16) 


s 


and 
if a* hp, 6,8 8, sin 9 a9 =— { ah&& (Q2.) sin@dé6; (17) 
0 8 


where Q, is the elastic potential energy per unit shell volume. Since 
motion occurs in two directions, there are two equations of motion. 
It is necessary to find expressions for R, and © so that both equations 
of motion are simultaneously satisfied. 

The radial equation, (16), is treated first. A value is specified 
for R, by equations (1), (8), (9), and (10). Combining equations 
(1) and (9), and integrating with respect to time, gives for the P, 
term of R,: 


n 


7A 
R,=«a ( —-—— qt P,,(cos 0) ef", (18) 


To find Q, it is necessary to know p,. Using (8), (4) and (12) gives 
us 
J On : 
p.— or (A, a” pi + Bn @™** po) Pn (cos 6) €7%"*. 


This can be simplified, with the aid of (10) and (18), and we get 


R, 

Ds] ps Wn? aC-—, (19) 

a 

where 
n 
Pian Po 
n+1 
c= (20) 
2 n Ps 


Likewise, replacing (12) by (13), gives us 
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R, 
Ds—J ps Wn” 03D, (21) 
a 
where 
n 
(ies Rupees 
f (Merge Ss 
D=—————_-. (22) 
N Ps 


Hooke’s law, applied in the radial direction, is used to eliminate 
e, from Q,. The generalized form of Hooke’s law applied to this 
problem gives us the following 
—p,=A(Q, + Ge + es) +2Ne1, 
where 4 is Lamé’s constant. For a shell we find 


Ata Dh 
Oy =—— (eat + €a? + @s*) + A (Cres + Cas + e061) 


ae z ie (€2? Jd} ( 1) 
exe + €2” + @3”) + C23 5 
2(A+ 2) H( 2 : ‘ a " 


where y = (1 + «)/(1—oc) and o=A/2(A + uw). Therefore 


Ds 
§ Q,=—— 8p, + ul (y—1)e + +1 6 @ 
Roe LL (y ) (y ) ps] 6 ee 


+ ul(y—Les+ (y + 1)e] dee. Ca 


Equations (14), (15), and (23) are used to eliminate e., e;, and Q, 
from (16). This gives us the radial equation of motion: 


Hs 2 bore t 
Qg— = =r aa 3 CO 6 ? 2 
a ‘ dé | ce 
where 
s2( 140; Anse ope eae 
a= sated Crmeet So Gece = 25 
h A+2u sh 
and 
9 Che ap, 
EI ae RT (26) 
bu 


: Equation (23) can also be used to expand the tangential equa- 
tion of motion, (17), in terms of e, and e,. After a partial integra- 
tion with respect to 6 of the term involving 0e./00 , and after elimi- 
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nating é, and e; by means of (14) and (15), the tangential equation 
of motion becomes: 


: d R; 
rast a 


d 1 d 
=— (2 +2)0,— (7 +) | ——="(@- Sin 6). 
dé sin 6 dé 


L 


(27) 


The two equations of motion of the shell, (24) and (27), and 
the expression for R,, (18), derived from the liquid motion, can all 
be satisfied simultaneously if, and only if, 


dR, 
O7—A ; (28) 
dé 
a=—2yn(n+1)A; 
and (29) 


2y=>([— (a+ 2) +n(n +1) (y'+1)]A; 


where A is a constant. Eliminating A from the last two equations 
gives us the relationship which specifies the characteristic frequen- 
cies: 


a[—(#+2) +n(n+1)(y+1)J=—4yn(n+ 1). (30) 


Finally, substituting the value of a from equation (25) gives us the 
equations for these characteristic frequencies in terms of x: 


bu 
2 u 


bu 
A+2u 


a 
2 —— 3 2 D— 2 
C? a8 [nme en 2 baer C+1+ a 


+[(1+0F) {amen +n —2| (31) 


tay femay ttt 1) — 2) 


The limiting case C = D = 0, reduces this last equation to the 
equation for the resonant frequencies in air. As a/h is large, the en- 
tire character of the roots has been changed. The cube term intro- 
duces an extra root. To arrive at numerical results it is assumed 
that « = 0.25 and therefore 4 = u. For simplicity, the choice 
po = pi = ps = 1 gm./em.? is also made. The roots of equation (31), 
for n = 2, then are: 


104 RESONANCE OF BIOLOGICAL CELLS 


a/h=10 a/h=5 
alz = 0.22 atz — 0.385 
ply = 14.5 plz = 14.9 
oly = 4.6 X 10° oly = 2.5 X 10° 
Ay = 0.24 alg = 0.24 
pAg = —6.3 pAs — —3.6 
cho —=—7X10* chs = —1.3 X 10°. 
The large value of ,Az for »%2 makes |h ,| © || , and in phase with 


R,. [This situation might be improved by altering equation (8)]. 


For the other two roots the assumption that |h e,| << |R| is valid. 
For both of the self-consistent roots of equation (81) the polar 
diameter expands as the equator contracts and vice versa. The tan- 
gential motion is at right angles to the equator. As the polar diam- 
eter contracts, the shell moves away from the pole in the solution 
represented by the first root, and toward the pole in the other one. 
Equation (26) is used to compute values of uw from the char- 
acteristic values of x plus the data found in the table. For P. cau- 
datum we find 
alo = 1.2 X 10? dynes/em.? for .%23 
olla = 0.8 dynes/cm.? for .%. 


Thus only the first root gives a reasonable value of pu. 

Conclusions. Two different cell models both lead to resonant 
frequencies for cells the size of paramecium within the observed fre- 
quency range. These different modes are illustrated in the figure. 
The different models and modes could be distinguished, in theory, by 
their harmonics. The experimental data are not at present suffi- 
ciently precise to locate these harmonics. 

No mention has been made of the actual cause of cell destruc- 


tion. It might be due to simple physical rupture, or to heat generated 
at the antinodes. 


SYMBOLS USED 


r,0,~= spherical polar coordinates (radius, colatitude and longitude 
respectively) 
R,®,¥= displacements from equilibrium position 
v= vector velocity of the liquid 
t= time 
f(t) = an arbitrary function of time 


“c= A dot over a symbol indicates partial derivative with respect 
to time 

Q = potential energy per unit mass 

® — velocity potential 

T — surface, or interfacial, tension 


EUGENE ACKERMAN 105 


p= pressure 
p= density 
P,,(cos 6) =the nth Legendre’s polynomial of variable cos 6 
== 27 times the frequency 
a= radius of the sphere 
subscript 7—= inside of sphere 
subscript 0 = outside of sphere 
subscript s = spherical shell 
é, = dilitational strain in r direction 
é, = dilitational strain in @ direction 
e,— dilitational strain in ¥ direction 
\== Lamé’s constant (A — wave length in equation on surface 
tension waves from Page) 
“= coefficient of rigidity 
o = Poisson’s ratio 
y= (1 + ¢)/(1 — 2) 
= 0? a? p,/p 
j= (1)? 
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FLUID MOTION 
NON-RIGID MEMBRANE ~ T mode 


Nodes at 0=55°, 125° 
Motion Radial Only 


Xp mode 
X, mode 


RIGID SHELL RADIAL 
and 
TANGENTIAL MOTION 


RIGID SHELL TANGENTIAL MOTION ONLY 


5 Radial Nodes at 0=55°and 125° 
fluid mode 
pamodeamaNe - he 4 Tangential Nodes ot O=0% 90° and 180° 
Nodes at 0=0°,90% and 180 Xq mode 


(e) 


RIGID SHELL RADIAL and TANGENTIAL MOTION 
Nodes Same As X, mode 


Xp mode 


Figure 1. The motion of the liquids and cell membrane or cell cortex. These 
motions are described in detail in the text. 
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The weak connectivity y of a random net is defined and computed 
by an approximation method as a function of a, the axone density. It 
is shown that y rises rapidly with a, attaining 0.8 of its asymptotic 
value (unity) for a = 2, where the number of neurons in the net is 
arbitrarily large. The significance of this parameter is interpreted also 
in terms of the maximum expected spread of an epidemic under certain 
conditions. 


Numerous problems in various branches of mathematical biology 
lead to the consideration of certain structures which we shall call 
“random nets.” Consider an aggregate of points, from each of which 
issues some number of outwardly directed lines (axones). Each 
axone terminates upon some point of the aggregate, and the prob- 
ability that an axone from one point terminates on another point is 
the same for every pair of points in the aggregate. The resulting 
configuration constitutes a random net. 

The existence of a path in a random net from a point A to 
a point B implies the possibility of tracing directed lines from A 
through any number of intermediate points, on which these lines 
terminate, to B. 

We shall say that B is t axones removed from A, if ¢ is the 
smallest number of axones contained in any of the paths from A to 
B. Point A itself is zero axones removed from A. All the other 
points upon which the axones of A terminate are one axone removed. 
The points upon which the axones from these latter points terminate, 
and which are not one or zero axones removed, are two axones re- 
moved, etc. 

The notion of a random net may be generalized, if it is not 
assumed that the probability of direct connection between every pair 
of points in the net is the same. In that case it is necessary to define 
this probability for every pair of points. This can be done, for ex- 
ample, in terms of the distance between them or in some other way. 
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If the connections are not equiprobable, we shall speak of a net with 
a bias. 

The following examples illustrate problems in which the con- 
cept of a net, defined by the probability of the connections among 
its points, seems useful. 


1. A problem in the theory of neural nets. Suppose the points 
of a net are neurons. What is the probability that there exists a path 
between an arbitrary pair of neurons in the net? If the net has bias, 
what is the probability that there exists a path between a specified 
pair? In particular, what is the probability that a neuron is a mem- 
ber of a cycle (i.e., there exists a path from the neuron to itself 
through any positive number of internuncials) ? Or, one may ask, 
what is the probability that there exists a path from a given neuron 
to every other neuron in the net? 


2. A problem in the theory of epidemics. Suppose a number 
of individuals in a closed population contract a contagious disease, 
which lasts a finite time and then either kills them or makes them 
immune. If the probability of transmission is defined for each pair 
of individuals, what is the expected number of individuals which will 
contract the disease at a specified time? In particular, what is the 


in Apa 1. The probability tree for the number of ancestors of a single 
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expected number of individuals which will eventually (after an in- 
finite time) contract the disease? Or else, what is the probability 
that the entire population will succumb? Note that if the probability 
of transmission is the same for each pair of individuals, we are deal- 
ing with a random net. 


3. A problem in mathematical genetics. Given the probability 
of mating between each pair of individuals in a population (as a 
function of their distance, or kinship, or the like), what is the ex- 
pected number of ancestors of a given order for each individual? 
Clearly, the less the expected number of ancestors, the greater the 
genetic homogeneity of the population. 


Each of these problems can be formalized by constructing a 
“probability tree.” As an example, a tree for the genetic problem 
is illustrated in Figure 1. 


We note that the tree consists of ‘‘nodes” connected by lines. 
The nodes can be designated by “first order,” ‘second order,” etc., 
depending on their distance from the “root.”” The number at the node 
indicates a possible number of ancestors of a given order. The lines 
connecting the nodes are labeled with the corresponding probabilities. 
Thus p,(2) = 1, since it is certain that an individual has exactly 
two ancestors of the first order (parents). However, the parents 
may have been siblings or half-siblings. Therefore it is possible that 
the number of grandparents is 2, 3, or 4. The corresponding prob- 
abilities are p.(2), p.(3), and p.(4). The probability of having a 
certain number of great-grandparents depends on how many grand- 
parents one has had. Consequently, those probabilities must be des- 
ignated by p;(i,7) wherei=2,----4and7—2,----8. In general, 
the probability of having a certain number of ancestors of order k 
will depend on how many ancestors of each of the smaller orders one 
has had. If, however, we simplify the problem by supposing that the 
probability of having a certain number of ancestors of the kth order 
depends only on how many ancestors of the (k — 1) th order one has,- 
then the probability that an individual has exactly n ancestors of 
the mth order will be given by 


am 8 4 
Pe (n) => sects Dy >, D2 (2,1) ps (1,7) Ds (3,k) a Pin (7,n) . (1) 
T=2 J=2 i=2 
The expected number of ancestors of the mth order will then be 


E(m) = nP(n). 29 
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Clearly, a similar tree can be constructed for the neural net 
problem. Here the numbers at the nodes of the kth order would des- 
ignate the possible number of neurons & axones removed from a given 
neuron. The p’s would designate the corresponding transition prob- 
abilities from a certain number of neurons (& — 1) axones removed 
to a certain number & axones removed, etc. If N is the number of 
neurons in the aggregate, clearly, a neuron B is at most N axones 
removed from a neuron A , or else there exists no path from A to B. 
Hence HE (N) represents the expected number of neurons in the aggre- 
gate to which there exist paths from an arbitrary neuron, if the 
neurons are not in any way distinguished from each other. This 
expected number we shall call the weak connectivity of a random net 
and will designate it by ». 

The contagion problem could be formulated in similar terms. 
Here weak connectivity would represent the expected number of in- 
dividuals which will contract the disease eventually. If we define I, 
the strong connectivity as the probability that from an arbitrary 
point in a random net there exist paths to every other point, then I’ 
will represent the probability that the entire population will suc- 
cumb in the epidemic described above. In this case, the number of 
“axones” represents the number of individuals infected by a carrier 
before he recovers or dies. 


The weak connectivity of a random net. We shall compute the 
weak connectivity of a neural net in terms of certain approximations 
whose justification will be given in subsequent papers. It will be 
assumed that: 


1. The number of axones per neuron a is constant throughout 
the net. This constant (the axone density) need not be an integer, 


since it may equally well be taken as the average number of axones 
per neuron. 


2. Connections are equiprobable, i.e., an axone synapses upon 
one or another neuron in the aggregate with equal probability. 


A. Shimbel (1950) has formulated the problem in terms of the 
following differential-difference equation 


dx/dt = [N — «(t)] [a (t) —a(t—7)]. (3) 


Here w(t) is a function related to the expected number of neurons ¢ 
axones removed from an arbitrary neuron, and - is related to the 
axone density. Then the problem of finding y is equivalent to the 
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problem of finding x(0«). A somewhat generalized form of equation 
(3) is given also by M. Puma (1939). The solution of the equation 
is, however, not given. 

An approximate expression for y where N is large was derived 
by one of the authors (Rapoport, 1948) where the number of axones 
per neuron is exactly one. This case will be generalized here to a 
axones per neuron, which are supposed constant through out the 
aggregate. 


The axone-tracing procedure. Let us start with an arbitrarily 
selected neuron A and consider the set of all neurons removed by 
not more than ¢ axones from A. Let x be the expected number of 
these neurons. Then evidently x = x(N, a, t) depends on the total 
number of neurons in the net, on the axone density, and on t. More- 
over, the weak connectivity of the net can be expressed as 


y(N ,a) =a(N,a,N)/N. (4) 


Since N and a are fixed, we shall refer to the expected number of 
points removed from A by not more than ¢ axones by x(t). Note 
that t is a positive integer. 

We seek a recursion formula for x(t) which will give us an ap- 
proximate determination of that function. To give a rigorous treat- 
ment of the problem, one would need to deal with distribution func- 
tions instead of expected values. For example, p(i,t), denoting the 
probability that there are exactly 71 neurons not more than ¢ axones 
removed from A, would determine the distribution for ¢. Succes- 
sive distributions (for t + 1, etc.) would then depend on previous 
distributions, instead of merely upon the first moments of these dis- 
tributions (expected values). The “probability tree’ method does 
take these relations into account. An “exact’’ approach to the prob- 
lem will be given in a subsequent paper. Meanwhile, however, 
we shall develop an approximation method in which it will be as- 
sumed that the expected value x(t) depends only upon previous ex- 
pected values, and, of course, upon the parameters of the net. 


The recursion formula. We now seek an expression for «(¢ + 1) 
— a(t). This is evidently the expected number of neurons exactly 
(¢ + 1) axones removed from A. We shall make use of the follow- 
ing formula, which may be readily verified. Let s marbles be placed 
independently and at random into N boxes. Then the expected num- 
ber of boxes occupied by one or more marbles will be given by 


N[1—(1—1/N)*]. (5) 
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In our axone-tracing procedure there are a[u(t) — x(t — 1)] 
axones of the newly contacted neurons to be traced on each step. 
Then the total number of neurons contacted on the (¢ + 1)th tracing 
will be, according to formula (5), 


N[1i— (1—1/N) s#-#@-1] , (6) 


But of these neurons the fraction x(t) /N has already been contacted. 
Hence the expected number of newly contacted neurons will be given 
by 


x(t +1) —a(t) =[N—a(t)][1— (1—1/N)ee-2e-n7],_— (7) 
which is our desired recursion formula. 


Determination of y. Let us set 


y(t) =N—<2(t). (3) 
Then equation (7) may be written as 
ub + 1) = y(t) A— 1/N) aera (9) 
or 
y(t +1) (1—1/N) 9 =y(t)(1—1/N) 4, © (10) 
Hence 
y(t +1)(1—1/N)*™ =constant=K . (11) 
We proceed to evaluate K. We have 
y(¢+1) =K(1—1/N), (12) 


But y(t) represents the expected number of uncontacted points in 
the tth step. Since before the tracing began one point constituted 
the set of contacted points, therefore we have 


y(0) =N—1, (13) 
and using formula (5), 
y(1) = (N—1)“" N-. (14) 
Letting ¢ = 0 in (12), we obtain 
K=N-* (N— 1), (15) 


Furthermore, since y(1) = y(0) and (1 — 1/N)~* > 1, we have 
y(2) S y(1), ete. so that y(t) is a non-increasing function of t 
(this is also intuitively evident from the definition of y). Since y=0 
for allt, y(t) must approach a limit as t grows without bound. Hence 
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Lim y(¢ + 1) =Limy(¢) = Y. (16) 
t>00 too 


Note that y = ~(N) may also be considered as Lim x(t) /N. This is 
t>0 


So since contacting no new neurons on any tracing implies that no 
new neurons will be contacted on any subsequent tracings. If we 
continue to carry out tracings “symbolically,” it is evident that at 
some tracing not greater than the Nth no new neurons will be con- 
tacted, and all subsequent tracings will be “dummy” tracings. 

Using equations (12) and (15), we see that Y satisfies the tran- 
scendental equation 


Soa (Ne) eNO. (17) 
For large N , this can be approximated by 

eee eraee Y/N 1) \ 1: (18) 
Hence, for large N , 

Y/N ~ Exp {a(Y/N —1)}. (19) 


But y= 2(0)/N=1—Y/N. Substituting this value into (19), 
we obtain the transcendental equation which defines y implicitly as 
a function of a , namely, 

y=1—e”™ (20) 


We note that for y = 0, every a is a solution of (20). If y #0, 
then equation (20) can be solved explicitly for a giving 


— log (t—) (21) 


a= 


Figure 2. Weak connectivity as a function of axone density. 
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The right side of (21) is analytic in every neighborhood of the 
origin and tends to unity as y approaches 0. Expanding that func- 
tion in powers of y, we have 


a=1+y/2+ 77/8 ----, (22) 


which allows us to plot a against y (cf. Fig. 2). This graph consists 
of two branches, namely, the entire a-axis and the function (21). 
Negative values of y, being physically meaningless, must be dis- 
carded. Thus in the region 0 S a = 1, we have y = 0, as is intuitively 
evident. We must show, however, that for a > 1, y follows the non- 
zero branch of the graph, otherwise we get the unlikely result that 
for sufficiently large N the fraction of individuals eventually infected 
in an epidemic will be negligible, regardless of the number of indi- 
viduals infected by each carrier of the disease. Actually, the solu- 
tion y = 0 is extraneous for a > 1 and appears in our equation be- 
cause we have let N increase without bound before determining the 
relation between a and y. In any physical situation N is finite. Hence 
a physically meaningful procedure is to determine y as a function 
of a and N and then allow N to increase without bound. Such a func- 
tion is given by equation (17). Proceeding from that equation we 
obtain 


Y/(N—1)=(1—1/N)*?™), (23) 

log Y —log(N —1) =a(N— Y) log(1—1/N), (24) 
= log Y —log(N — 1) 

"~ (N— Y) [log(N—1) —log(N)] 


Let us write Y = N — 6(N) = N[1—¢(N)/N]. Then equa- 
tion (25) may be written as 


(25) 


log N —log(N —1) + log[1— ¢(N)/N] 
=a¢(N) [log (N—1) —log N]. a 


Since ¢$(N) < N for all N, we may expand the last term of the left 
side of (26) and obtain 


log N — log(N — 1) —¢(N)/N —4[¢6(N) /N]? 
—1/3[¢(N)/N]*----=a4(N) [log(N —1) —log N]. 


We now expand log(N — 1) — log N which appears in the right side 
of (27) and after rearrangements obtain 


(27) 
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log N — log (N — 1) 
N. = ee 
=2O 1a AM 8, WoT a 


eee 2 
2N 3N2 | (28) 


< a if 
= | Sat 900) 0) 


Now if a is fixed and greater than unity, the limit of $(N)/N can- 
not be zero as N increases without bound, because otherwise for N 
sufficiently large the right side of (28) becomes negative, while the 
left side is always positive, a contradiction of inequality (28). There- 
fore, the limit of Y/N, as N increases without bound, cannot be 
unity fora > 1. But this means that y # 0 if a > 1. Hence, for 
a > 1, the non-zero branch of our curve is the only meaningful one. 

An examination of the meaningful part of the graph of equa- 
tion (20) shows that as long as the axone density does not exceed 
one axone per neuron, y = 0, ie., for very large N , the number of 
neurons to which there exist paths from an arbitrary neuron is neg- 
ligible compared with the total number of neurons in the net. On 
the other hand, as the axone density increases from unity, y increases 
rather rapidly, starting with slope 2. Already for a = 2, y reaches 
about 0.8 of its asymptotic value (unity) and is within a fraction 
of one per cent of unity for quite moderate a (say > 6). This means 
that no matter how large the net is, it is practically certain that there 
will exist a path between two neurons picked at random, provided 
only the axone density is a few times greater than unity. The in- 
terpretation in terms of an epidemic with equiprobable contacts is 
entirely analogous. 


The case a= 1. This case was treated by one of the authors 
(Rapoport, 1948) by a different method. It was shown that for large 
N, the probability that a neuron was member of a cycle was given 
by V/x/2N. This gives the probability of the existence of a path 
from a neuron over any number of internuncials greater than one to 
itself. But under the assumption of equiprobable connections, this 
may well represent the probability of the existence of a path from the 
given neuron to any other neuron in the net. Therefore we should 
have for large N, in the caseea=1, 


y~ V/n/2N . (29) 


For N = o, y reduces to zero, as it should according to equation 
(20). We shall, however, examine the asymptotic behavior of y for 
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large N deduced from our approximate method, in order to compare 
it with the asymptotic behavior (29) deduced from an exact treat- 
ment of the special case. Dividing both sides of (17) by N, we may 
write fora=1 


Y/N=[(N—1)/N]*™™, (30) 
whence, since Y/N =1—y, 


day GN Na 
(31) 
= Exp{In(1—1/N) + Nyin(t—1/N)}. 
We let z = N~ and examine the behavior of y for small values 
of z. Expanding the right side of (31) by power series and retain- 
ing only terms of the second order (note that z and y vanish to- 
gether), we obtain 


Ll—y=14+[—-2—822-:-- | + [yy — 2/2 —- *] 


(32) 
PRA chee ye 2 FV is aioe 
Hence, 
O=—2 4+ y?/24+ y2/2+----, (33) 
Differentiating with respect to y, we get 
dz/dy=y + y/2.dz/dy + 2/2 +--+, (34) 
dz/dy ~ (y + 2/2)/(1— y/2). (35) 
Therefore dz/dy vanishes at z = 0, y = 0. Differentiating once 
again with respect to y, we obtain 
d?z aa (36) 
dy? 2=0 : 
y=0 


Hence the power series representing z as a function of y begins 
as follows: 


Z= 92/2 eee, (37) 
Thus 

y? ~ 22=2/N , (38) 

y~ V2/N =141 VN. (39) 


The “exact” result as expressed by (22) gives 
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y~1.2/V/N. 


Thus the approximate method applied to the case a = 1 implies an 
asymptotic behavior of y for large N which does not depart too 
sharply from that deduced by the exact method. The limiting value 
for y is zero in both cases. The question of how well the limiting 
values of y are approached by the approximate method for a > 1 
remains open. 

This investigation is part of the work done under Contract No. 
AF 19(122)-161 between the U. S. Air Force Cambridge Research 
Laboratories and The University of Chicago. 
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This paper deals with diffusion into a medium of finite thickness 
through a flat structure which can be considered either as a slice of 
tissue or as a membrane. Formulae are given to determine the diffusion 
coefficients in both the flat structure and the adjacent medium from the 
knowledge of the amount of substance penetrating the medium. The 
meaning of the formula: (permeability coefficient) — (diffusion co- 
efficient) /(membrane thickness) and the experimentally observed vari- 
ability of the permeability coefficient in the non-steady state are in- 
terpreted on the basis of the mathematical theory of diffusion. 


1. The background and the objectives of the paper. There are 
essentially two methods for the determination of the permeability 
coefficient. One whose theory was given by D. H. Andrews and J. 
Johnston (1924) consists of immersing the membrane or the tissue 
slice in the substance for which the permeability of the tissue is be- 
ing studied and of measuring the gain in weight of the slice in rela- 
tion to the duration of the immersion. This method has been used 
both by chemists and physiologists, particularly in the case of liquids 
where the gain in weight is more easily determinable than in the case 
of gases. For the latter the method developed by H. A. Daynes 
(1920) has been preferred. It consists of letting the gas diffuse 
through the membrane rather than into it as in the previously men- 
tioned method. It can be diagrammatically represented as shown in 
Figure 1. The tissue slice or membrane is bounded by the planes 
2—=—aand x=0. The gas is kept at a constant concentration C 
in the region « < —a and diffuses through the membrane into a con- 
tainer bounded by the planes x = 0, x = b. The initial concentra- 
tion of the gas in the slice and in the container is assumed to be zero 
and the permeability coefficient of the slice is determined from the 
relationship between the amount of substance accumulating in the 
container and the duration of the diffusion. The theories of both 
methods consider the diffusion to be a one-dimensional process. An- 
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FIGURE 1 


drews and Johnston assume the concentrations to be equal on both 
faces of the slice and constant throughout the experiment, an as- 
sumption which is well satisfied if the thickness of the slice is small 
with respect to the mass of the surrounding liquid. Daynes assumes 
different concentrations on the two faces of the slice but constant 
in time, that is, unchanged by the amount of substance diffusing 
into the container. This latter condition requires either a limitation 
of the experiment to negligible changes of concentration in the con- 
tainer or the provision of a special arrangement which maintains 
constant the concentration in it. For instance, the gas diffusing into 
the container may be removed or its partial pressure may be kept 
constant in any other manner. In the usual applications of this 
method the choice is made to fulfill the assumptions of the theory by 
limiting the measurements to small changes of concentration in the 
container. : 

It is the objective of the present paper to show how this limita- 
tion can be avoided, so that measurements covering large ranges of 
concentration can be used. An experimental arrangement which cor- 
responds to the conditions of the present theory was designed by H. 
Miiller (1941). The present paper gives methods for a determination 
not only of the membrane permeability, but also for a simultaneous 
determination of the diffusion coefficient in an adjacent medium. 
This is more of physiological than of purely chemical interest, since 
artificial membranes are usually obtainable in isolated form, but not 
so the living ones. The methods given include one which does not re- 
quire any determination of concentration in, or weight of, the mem- 
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brane, a significant detail from a practical point of view. 

2. An exact solution of the problem. Although we refer to the 
media bounded by the pairs of planes ¢ = — a; x = 0 and xz = 0; 
x = b as membrane and container respectively, various interpreta- 
tions are possible. For instance, the two media can represent two 
adjacent layers of a tissue or a tissue layer and an adjacent part of 
animal body. 


Let C.(z, t) and C,(x, t) be the concentrations of the substance 
considered in the membrane and in the container respectively. Let 
D, and D, be the corresponding diffusion coefficients. If vacuum ex- 
ists initially in the container and the diffusing substance is a gas, D, 
is its self-diffusion coefficient. The following equations rule the dif- 
fusion of the substance from the medium where wz < — a into the two 
media between — a and b under the assumptions that the substance 
was not present initially between — a and b, that its concentration 
at « = — a is always kept constant and equal to C and that no sur- 
face phenomena exist on the membrane which could disturb the con- 
tinuity of concentration at the boundaries: 


0 Cx 0? Ci 6G. 0C, 
ee — D; =), 

at 0 ax OD |e oe (o 

oC 

eee Oe 0,. CH — a4!) —C. 
Ox a=b 


C.(0, t) =C;(0,¢), 


where k =a, b. We solve this system of equations by the method of 
the Laplace transformation. Using the following symbolism: 


f(s) =L{ F(t)}= fret Fat, F(t) =| f(s) 


and putting c,.(“, s) = L{C,(z, t)}, we obtain: 


SC; = D; 0? Cy /0 Xx? , D,0 Cy/0 £20 = Dy 04/0 %] e-0 » 
0¢,/0%]e.=0, G(—a,8) =C/s,. ¢,(0,$) =%(0,8). 


Putting, as before, k =a, b and 


ES eV S/ Lk > Oe — ky/s/Dr y DSV0/Ds; 
6=D coth o-cosho, + sinha , 
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the solution of the above system of equations is: 


€, = C(D coth s - cosh &, — sinh &,)/(6 $) 
¢, = DC (coth o; - cosh & — sinh &)/(6 s). 


Consider a constant cross section having an area A. The amounts 
of substance which at t = o are present in the membrane and in the 
container are ACa and ACb respectively. Let Q.(t) and Q,(t) be 
the fractions of these amounts which are present at any moment ¢. 
Let us introduce the notations 


Qa(s) =L{Qa(t)}; a(s) =L{Qz(t)}. 


Then we have 


(8) =L(Q.(t)} = (a0)> f caler, t)de 


- -4 


= (D,/s)1(D coth o- sinh o, + cosh o,—1)/(a6ds), 
qo (8) = L{Qz(t)} = (bC) (e oy (a , t) dx = (D,/s)¥2/(b 8s). 


Since the explicit expressions of Q,(t) and Q,(t) are very compli- 
cated and since Q,(t) is of less practical interest than Q,(t), as far 
as the objectives of this paper are concerned, we limit ourselves now 
to the consideration of Q,(t) only. The mathematical procedure used 
to obtain Q,(t) can also be used to obtain Q,(t). 


If we set A= (1— D)/(1 + D),c=o, + % and 
R= ([1—A exp(2 o,) — A exp(2 o)] exp(— 2c), 


we have 


qo(s) = 4[b(1 + D)s(1'+ R)]> VDyz/s exp(—) sinh o,. 


If s is sufficiently large R is numerically smaller than unity. Con- 
sequently the binomial expansion can be applied to 1/(1 + R). Using 
the formula 


(1+a+ b)*=3 nlar itt ae Yulyt] , 


by 
where the sum = is taken for all non-negative integers u and » such 
that 0 < w+» < 1, we obtain: 
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SS dette 


eo 


= 2 D Bayy- exp[—2(n— w)o.—2(n—v) op), 


n=0 pL,v 
where 
Bnpjv = (—)™ 2! AMY/ (mn — w—v) Yuly!] 
Also, putting 
Anpy = (2n— 2u +1) aD? +2(n—») dD, 
Kayv = Any + 2bD,, 


we obtain 
a» (s) = 2[b(1 + D) 8] (D,/s)** SS Bo ys [exp (Hane V3) 
— exp (—Kiyy Vs)]. 


From the above we have, using a known formula of the Laplace trans- 
formation (Churchill, 1944, p. 299, formula 85), 


Q(t) =2[6(1 + D)]7* V Did 2 Bowel Bi (05 nn) 
— E(t; Knyy)], 
where, as previously, the sum 2 is taken for all non-negative RES 
# and y such that 0 < w+ mand 
«E(t; K) =K[(r Va) exp(— #2) —erfe(r)] 
with 7 = K/(2 vt) and 


2 +00 
erfe (7) =—— ip exp (—a?) da. 
Vit es 


Thus E can be calculated by means of two well-known functions and 
also by means of a single function tabulated by D. R. Hartree (1935) : 


ierfe(r) = i; erfe(x) da. 
ao 
In fact, H(t; K) = (K/r)ierfc(z7). 
From the MacLaurin series: 
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erfe(s) =1— (2/Va)r+-: (1) 
we obtain the following for small 7: 
ierfe(r) = (1/Va) —1 + (P/V a) tee (2) 
From the above, for large t, we obtain 
E(t; K) = 2Vt/n—K + [K?/(2Vat)]. 
CPUS state aco 


fee) 


Q, (t) =? 2[b (1 sr D) 12 Waal = 2 Bap, (Kayp,v Saale Anu,v) . 


Using the expressions of Bnyv, Kny,v, and Any» it can be checked 
that the right-hand side of this last relation is unity, as it should 
be since Q,(0) =1. In addition E(0; K) = 0, so that the condi- 
tion Q,(0) = 0 checks also. 


3. Approximate methods. Since the calculation of Q,(t) by 
the formula previously derived is, in general, very laborious we give 
approximate expressions of this quantity as well as approximate 
methods for the determination of D, and D, from experimental data 


on Q,(t). 


(i) Q.(t) for large values of t. We derive a formula for 
Q;(t) valid for t > o by applying the following Tauberian theorem 
of the Laplace transformation (Doetsch, 1937, p. 208). If F(t) > 0 
fort >0,if L{F(t)} =f(s) is convergent for s > 0, if 


lim [A (ux) /A(a4)] =1 


for any wu > 0 and f(s) ~ A(1/s) for s > 0, then 
D> {s4f (s)} —1 () for f= 


where ~ is a symbol of asymptotic equality in the sense that the ratio 
of its two sides tends to unity as the independent variable tends to 
its limit. Take s*f(s) = q(s) and consequently f(s) = sq)(s). The 
applicability of the Laplace transformation to the system of equa- 
tions treated here requires the existence of I-{sq.(s)}. In fact, 
the existence of L{dC,/dt} implies the existence of L{dQ,/dt} and 
the latter equals sq,(s) . However, since we have not even discussed 
the applicability of the Laplace transformation to our problem, as- 
suming it as justifiable by physical reasons, it is worthwhile to show 
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at least that L-*{sq,(s) } exists, a detail of importance for the applica- 
bility of the Tauberian theorem just stated. 
Since 
8qo(s) =b* VD,/s (sinh o, sinh «, + D cosh o, cosh «)- sinh a; 


and the term in parentheses of this expression considered as a func- 
tion of s has all real and negative roots (Churchill, 1936), sq,(s) is 
an analytic function in any region of the complex s-plane in which 
the real part of s is > 0. The fact that s enters in the above ex- 


pression through \/s does not introduce any branch point, since 


(sinh o)/\/s, sinh o, sinh o, cosh o,, cosh o are all single-valued 
analytic functions of s. In addition to this since, with e any finite 
positive number, 


lim s?** q,(s) —=0, 
I{sq,(s)} exists for s > 0 (Churchill, 1944, p. 159) and equals 
dQ, (t) /dt. Consequently, in applying the Tauberian theorem stated 
above we may take F(t) = dQ,(t)/dt. This function satisfies the 
requirements of that theorem. In fact dQ,/dt > 0 because Q,(é), is 


by its physical meaning, a monotonically increasing function. As it 
has just been shown L{dQ,/dt} is convergent for s > 0. Also, put- 


ting 
en a Wee 
Sara 2 BD,’ 


as a a’b? bt 
aay ( Friis "6D.D,  45D\"" 
the following expansion is valid for small s > 0: 
Qp = 87/(1 + fis + B.S? +--+) = SZ, 
where Z =1— fs + (B12 — B2)s? + ----. Therefore: 
L{dQy/dt} = sq (s) =Z 
and in applying the Tauberian theorem we can put 
f(s) =sq(s) ~ A(1/s) 


with 4(s) = 1 — fis? + (f,? — B2)s”, because for this function 1, 
the condition 


lim [A (us) /A(s)] =1 
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is satisfied. Therefore, fort 7 o, 
L~{sf (s)} = Q,(¢) ~1— pit * + (Bi? — Ba) o- (3) 


This considered as an equation in Q, and 1/t represents an arc of a 
parabola. Fitting it to the experimental data the constants 6, and fz 
are obtained and from them 1/D, and 1/D, are calculated as solu- 
tions of a system of two algebraic equations, one of first and the other 
of second degree. Of course one should not expect to obtain through 
this procedure more than approximate values of the diffusion co- 
efficients D, and D,. However, these can be checked and corrected 
by substitution into the exact expression of Q,(t). The question 
may arise: How much significance should be attached to this appli- 
cation of the Tauberian theorem since from a purely mathematical 
viewpoint equation (3) implies only the statement that the limit 
of the ratio of Q,(t) to the right-hand side of equation (3) tends to 
unity as ¢ tends to «? Consequently, changing (6, and f, in equa- 
tion (3) into any other finite numbers would give an equation which 
would be asymptotically still correct, but which would imply differ- 
ent values of the diffusion coefficients D, and D,. The answer to this 
question is as follows: The great majority of applications of approxi- 
mate methods in mathematics is to some extent empirical and must 
be so since estimation of errors is, in general, impractical or impos- 
sible. The type of approximation involved in the above application 
of the Tauberian theorem is essentially the same as that inherent in 
the application of the MacLaurin series, which is the most popu- 
lar of all approximations. In fact, saying that the power series 
f(x) =1+ f'(0)x% + ---- can be approximated for a given small x 
by its first two terms is often equivalent to saying that 


lim f (2) /[1 + f (O)e]—=1. 


But this equation would be true even if f’(0) were changed into any 
other finite number. Nevertheless, the common linear approxima- 
tion of a function implies the use of f’(0) as the coefficient of x. The 
reliability of an approximate method is, in general, not proven a 
priori but tested through application. It will be shown in the next 
subsection and also in section 4 that completely different approxi- 
mate procedures lead to formulae which are practically identical with 
those obtained through the above Tauberian theorem. This speaks 
for the reliability of the latter. 

A procedure similar to the one applied for Q;(t) is possible also 
for Q.(t) if data on the amount of substance existing in the mem- 
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brane are available. One detail of this procedure deserves explicit 
mention. We have, for small s > 0, 

Qa(s) =s*(1—as+....), 


where a = (a/D.)[(a/8) + (b/2)] and through an application of 
the Tauberian theorem in a similar manner as before we obtain for 
Geo? 


Qu (t) ie at ot 


Thus for large t the diffusion into the membrane is independent of 
the diffusion coefficient D, in the container. The above formula gives 
a very simple relation for the determination of an approximate mag- 
nitude of the diffusion coefficient D, in the membrane or slice of tissue 
from experimental values on the amount Q,(t) of substance exist- 
ing in it. But the measurement of Q,(t) is often impossible or diffi- 
cult. 


(ii) Another approximation. It has been shown (Opatowski 
and Christiansen, 1950) that exp (—gr) with g equal to about 1.6 
is an approximation to erfc(r) in the range 0 to + o. Consequently 


ierfc(r) = g exp(—gr), 
E(t;K) ~ K(g 7)“ exp(—g r) with> = K/(2 V2). 


These relations are exact for 7 = o. An idea about their degree of 
approximation for small 7; is offered by the following comparison: 


ierfc(0) =1/\/a = 0.564, 
[og exp (—g 7) Jroo ~ 1/1.6 = 0.625. 


Let us now consider a function f(s) which admits the following 
expansion with k, > 0: 


f(s) =¢(Vs)/(s Vs) with $(Vs) a A, exp(—kn V8). 
Since L{E (t; K)} = (s Vs) exp(— K v/s) we have, under the as- 


sumption that the inverse Laplace transformation can be applied 
termwise to the above sum >, 


F(t) =L“{f(s)} é 
=> AE (t; kn) = (Vt/7) = An exp (—y kn/ Vt), 


where y = g/2 = 0.8. But this latter sum is obtainable from that 
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equal to ¢(\/s) through the substitution V/s > y/V/t. Consequently, 
F(t) =L*{o(Vs)/(s V8)} = (Vt/7) ¢(r/V4). 


It is seen that this approximate formula, whenever it is applicable, 
gives formally the same result as the Tauberian theorem previously 
indicated, except that ¢ is replaced by t/\/y ~ 1.12¢. Whereas the 
Tauberian theorem gave a formula valid for large ¢, the present 
method yields the following approximate result valid for any ¢: 


Vv Dit by p ay 
Q,(t) = —— D coth —  )cosh 


ay 
+ sinh a - 
(0s) 


A similar expression holds for Q,(¢). 

A problem of the diffusion theory in which a constant concen- 
tration is assigned at x = b (see Fig. 1) instead of the condition of 
zero flow at « = b has been discussed by Churchill (1936). It can be 
treated mathematically in the same way as the problem here consid- 
ered and would yield an extension of the immersion method of An- 
drews and Johnston for permeability determination to the case of a 
tissue slice consisting of two different layers. This method could be 
further extended to multi-layer tissue slices; some mathematical work 
usable toward this end is available (Jaeger, 1950). 


4. Approximation based on a different diffusion problem. We 
give in this section an approximate method for the determination of 
the diffusion coefficients D, , D, based on known formulae of the dif- 
fusion theory. These formulae deal with the case in which the diffu- 
sion through the membrane takes place into a semi-infinite medium, 
that is, when b = o. We have then (Carslaw and Jaeger, 1947, pp. 


261-264) : 
a (Q2n+1l)at+e2 2n+ 1l)a— 
Ca=C > A erie eric GP Sora 
see 2V Dit 2V/ Dat 


3 


2 (2n+ 1)a+ 2D - 
C,=C(1— A) > A* erfe ———_———__ 
n=0 2VD.t 
The amount of substance existing in the container at a time t ex- 
pressed as a fraction of the amount that would exist at t = o can 
be calculated approximately by the following formula: 
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Qs(t) = (bC)> f O@.pae. 


Since the concentration in the container is determined by free diffu- 
sion, the amount of substance diffusing during any time into a con- 
tainer of finite length is smaller than the amount that would diffuse 
if that container had an infinite length. The integration in the above 
formula can be carried out by means of the function ierfe already 
considered. However, we limit ourselves to an approximate expres- 
sion using relation (2). We obtain in this way for large t: 


E Db — 2a | 
2\/n Dat Va Dit 


Q(t) =(1— A) > A” 


n=0 


If we restrict ourselves to the case, most frequent in applications, in 
which D, < D,, we have 0 < D<1land0< A< 1. Therefore we 
can use the expansions: 


3 A"=1/(1— A), 3 (n+ 1) 4"=(1— A) (5) 


n=0 


and obtain from the previous expression of Q,(t), for large t, 


Q, = 1— (a/v), (6) 


1 /;avD; b 
Wp — Tae at ees 5 
Da 2D, 


where 


In a similar manner we obtain for large t 
0 
Qa(t) = (aC) i C.(a,t)dx =1— (a/Vt), (7) 


where @, = aVD;/(2D.\/x). Thus by fitting experimental data for 
Q,(t) and Q,(t) with the above asymptotic relations, the diffusion 
coefficients could be approximately calculated. The explicit formulae 
are: 


b? a Dz 
Ds De= oi 
An (204 — wr)? 2a mt 


If the diffusion coefficient D, in the container is known a priori, then 
the experimental determination of Q,(t) alone is sufficient to give 
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the diffusion coefficient D, in the membrane. In fact, the equation 
for w;, solved with respect to D., gives the latter in terms of known 
quantities. 

The same approximate formulae can be obtained also by apply- 
ing the following Tauberian theorem: If F(t) is > 0 and monoton- 
ically decreasing for t > 0, if L{F(t)} = f(s) is convergent for 
s > 0, and 


f(s) ~ Bs otorgaiee0, 
with 0 < y <1, and K a constant > 0, then 
F(t) ~ Kit] (y—1)! > for, £7 cs. 


G. Doetsch (1937, pp. 208-09) gives a proof of this theorem when 
F(t) is a monotonically increasing function. The proof for a mo- 
notonically decreasing function can be given in a similar manner. By 
putting coth o = 1 since b = o, Ga, c, and q, are now obtainable 
from the formula of section 2 (see also Carslaw and Jaeger, 1947, pp. 
261-64). Consequently expanding in series we obtain: 


Qa(s)~s2?(1—a@aVus) for s>0. 
Since P,(t) = 1 — Q,(¢) is by its physical meaning a monotonically 
decreasing function and 
Da(s) =L{P.(t)}~ a@aVa/s for s>0, 
we have, using the Tauberian theorem just stated, 
P(t) = of Cb for BES oe 
which leads exactly to equation (7) for Q.(t¢). If, instead of this 


Tauberian theorem, the one applied in section 3 were used a formula 


like (7) would be obtained with the difference that w.\/a would ap- 
pear therein instead of w,. Thus the approximate character of the 
result would be essentially the same. Similar conclusions are ob- 
tained for Q;(t¢) if we use 


b 
a (s) = (BC) f co(w , 8) da 


= (bs) VD,/s [1— exp(— 0) ]/(D cosh o, + sinh o,). 


5. Diffusion coefficient in the membrane and membrane per- 
meability. Some of the formulae derived in the previous section 
make it possible to discuss the relationship between the diffusion co- 
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efficient D, in the membrane and the permeability of the latter. We 
take as the mathematical expression of the permeability: 


h=— Dy(C— Cy) 8Cy/02| a0 


which equals the number of molecules leaving the membrane per unit 
area and unit time when the concentration difference on both sides 
of the membrane is unity. Actually the common definition of per- 
meability concerns the number of molecules that pass through the 
membrane, but it is difficult to express this mathematically since the 
instantaneous time rates of flow are, in general, different at both 
faces of the membrane. Substituting the expression (4) of C, into 
the above formula for h it is easily seen that h changes with t , and 
this has been experimentally observed even on artificial membranes 
of simple chemical composition (cf. e.g., Tolliday, Woods and Har- 
tung, 1949). However, it can be shown that h > D,/aast > o. 
In fact, when ¢ is very large 


i) C,/0 | no =~ C/\V/a Dyt 
and, using equations (4), (1) and (5), we find 


Cc Bite eiteeatnt (pus. 02 ) 
Ine = 6 | al tp 


From here, taking into account the expression of 4, we obtain 
h = D,/a, which is a known intuitive relationship. The relation of 
inverse proportionality between permeability and thickness has been 
experimentally observed on artificial membranes (Sager and Sucher, 
1939; Edwards and Pickering, 1920; Daynes, 1920; Carson, 1934). 
In the physiological field the observed relationships appear often to 
be more complicated (see e.g., Krogh, 1919) ; one reason for this may 
be the dependence of the diffusion coefficient on concentration. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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A formula is derived for the probability of hitting exactly k new 
targets with s shots, where it is supposed that m out of N targets have 
already been hit and that each shot results in a hit. 


Suppose one shoots at random at an area in which N targets 
are closely packed so that each shot is a hit at some target, where the 
probability of hitting a specific target is 1/N. Then it is easy to 
verify that the expected number of distinct targets to be hit with s 
shots will be 


BN, 8) N fie 1/N) 2]. (1) 


If m targets have been hit, the expected number of new targets 
to be hit with s shots will be 


(N jm) = (N —m) [1-41 =1/N)*] . (2) 


These results are useful in the theory of random nets (cf. Solo- 
monofi and Rapoport, 1951; Rapoport, 1951). However, it is also 
desirable to know the probability that exactly k new targets will be 
hit with s shots, where m targets had already been hit, for each 
k(0 = k = 8s). We shall show that this probability is given by 

an ! ke 
a ee (-1)/( * m+ k—a'. (3) 
(N—m—k)! Kk! N® j20 J 

As we fire our s shots consecutively, one of the following mu- 
tually exclusive events will occur: either the first new target will be 
hit on the first shot, or on the second, or on the third, ete. Whenever 
this event occurs, we shall hit exactly k + 1 new targets in all, pro- 
vided we hit exactly & targets with the remaining (s — 1) or (s — 2) 
or (s — 8) -:-- or & shots. These conditions are translatable into a 


recursion formula for the r’s , namely, 
133 
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m m 
Tra (S,m) = ge desi) tees a 


4 
m \2 Mm 8-k-1 (4) 
+( SY ne—simen t-( a rilesm+ 1) ‘ 
N N 


m 
On the right side of (4), the powers of ( eo represent the 


probabilities of the corresponding number of successive misses; 
(N — m)/N is the probability of the first hit; the arguments of % 
indicate how many shots are left and that in all m + 1 targets have 
been hit after the first new hit. 


For k = 1, the right side of (3) reduces to 


N—m 
Fi (8s 1) esas ee eet "padtlel (5) 


On the other hand, straight-forward computation of 7,(s, m) gives 


N—m[/m+1\* m mri\™ 
10 (8 fl) = eee | oe ne Bese: 
N ( N ES N | 
f1 N—-m 
ait 3) [aaa lore De mem De (6) 


N—m 
eee m4] as a [(m + 1)*§—™m']. 


Thus (3) holds for k = 1. We shall prove our result by an in- 
duction on k. To prove this induction we shall need the following 
lemma.* 


Lemma. If m and n are positive integers and m < n, then 


3 (5) iso 3(G)er=cyem. oD 


Proof. Let u(z) = (1 — e*)™. Then all the derivatives of u(x) of 
order less than n will involve only terms containing (1 — e”) as a 
factor, while the nth derivative 


wu™ (x) = (—1)” nle™ + terms containing (1 — e?). (8) 


*The proof of this lemma was communicated to the author by H. G. Landau. 
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Hence 
u™ (0) =0 form <n; u™ (0) = (—1)"n!. (9) 
On the other hand, 
> : 
u(x) meal 1 )er+ (9 ome + (a-1) 26%, (10) 
ny 
u® (¢) => ; Jay je, (11) 
j=0 


so that 


am (0) = 2 ( i) (—1)? jm: ym (0) =i ) (—1)/ q”. (12) 


Substituting (12) into (9), we obtain (7). 
Corollary. If r is any positive integer, and m < n, then 
Epp aes 
= (5) aid+ne=o. (18) 
f= 
Proof. Evidently it is sufficient to prove the result for r= 1, since 
in that case (13) holds for arbitrary r by induction. We have 


G+D9=3(7)m+ ig (14) 


Then (13) can be written as 


m1 n nN ; 
=(7)/ 5 (5) ar | +3(G) (15) 
~ 4=0 a j=0 Jj , j=0 J 
But the first term vanishes by our lemma, and the second term, being 
equal to (1 — 1)” vanishes also. Hence the corollary is proved. 
We now proceed with our induction on k& to prove (3), and as- 
sume that (3) holds for 7, with arbitrary arguments. 
Substituting expressions of the type shown on the right side of 
(3) into (4) we obtain the following expression for Ti: (8 , ™m) 


(N—m)! 


(N—m—1—k)!(k +1)! 
X {Lm +k + 1)84—k(m + bye tee tb (1) (m+ 1) 27] 


-s 


+m [(m +k +1)82—k(m + hye? +--+ (1) *(m + 1) 27] (16) 
+m? [(m +k +1)**—k(m + ke + +--+ (—1)¥(m + 1)**] 
- t ° e ry e 


+ me [(m + e+ 1)*—Kem + Be) eves + (1) (m+ 1)". 
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We now proceed to sum the terms within the braces of expres- 
sion (16) by columns. We observe that if the columns contained 
(s—1) terms each, ie., if the powers of m outside the brackets 
stretched all the way to m*", the columns would be of the form 

> ay? = (4° —y*)/(4—y). (17) 


a+b=s-1 


Therefore we add and subtract the missing terms so as to take ad- 
vantage of formula (17). This gives 


(N—m)! (m+ keel)? =a 


(N=m 1 =f) ea)! k+1 
. (m+ k)s—mt | k(e—1) (m+k—1)*—m 


—s 


2! kA 
(oD aes eae ERED é 
i (N—m—1—&k)!(k+)1)! 
x} Lm*(m +h +1)** + ms (m +k +1) (18) 


+.... me-*| — kL ms*(m aL k)* + ms-e (m + k)* ak oe me] 
k 
+( ) bm-*(m + 1)" + me (m + b— 1)? $e me] 


+ (—1)"* [m*(m + 1) + me (m +1)? +.--.m]!. 
Consider first the terms within the first brace. Factoring out 
(k + 1)-+, we may write them as 


N— 
( et Ne ([(m +k + 1)*—m'] — (k+1)[(m + k)*—m'] 


(k+1)k 
ear [(m + k—1)*—m'] Seen 


+ (—1)*(k +1) [(m+ 1)*§— mJ} 


(19) 


£ (ear) 4 5 (* 54) (m+ eam 
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But by a well known result on binomial coefficients, 


ere 
143 ("5 \p= ca. (20) 
Therefore the right side of (19) reduces to 
N—m el 
EH pee 1) ( ; J(m+e+1—3), (21) 
which is of the same form as (3). Our proof will be complete, if we 
show that the expression within the second brace of (18) vanishes. 


Let m + 1=~,r and note that each column but the last within 
the second brace of (18) is of the form 


k 


Em(5)oigen, (22) 


and this sum vanishes by our corollary above. The last column also 
vanishes in view of the argument following equation (15). This 
completes the induction. 

For the special case m = 0, k = 3, the right side of (3) should 
reduce to s! = k!, since an elementary calculation shows that the 
probability of hitting initially s targets with s shots is 

N! 
4 (CA eh oer rn ee (23) 
(N —s)! N° 

This result is indeed obtained if one substitutes the second equa- 
tion of (7) into (3) with m = 0. This is a check provided for the 
formula. 

H. G. Landau has also suggested an alternative preot of (8) 
based on an induction on s. 

It can be verified that 


s i epny = (NE) fl AN) (24) 


as should be the case, since the right side of (24) represents the ex- 
pected number of newly hit targets. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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